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Abstract. It is shown that, when the inflaton field modulates the gauge kinetic function 
of the gauge fields in supergravity realisations of inflation, the dynamic backreaction leads 
to a new inflationary attractor solution, in which the inflaton's variation suffers additional 
impedance. As a result, slow-roll inflation can naturally occur along directions of the scalar 
potential which would be too steep and curved to support it otherwise. This provides a generic 
solution to the infamous eta-problem of inflation in supergravity. Moreover, it is shown that, 
in the new inflationary attractor, the spectral index of the generated curvature perturbations 
is kept mildly red despite eta of order unity. The above findings are applied to a model of 
hybrid inflation in supergravity with a generic Kahler potential. The spectral index of the 
generated curvature perturbations is found to be 0.97 - 0.98, in excellent agreement with 
observations. The gauge field can play the role of the vector curvaton after inflation but 
observable statistical anisotropy requires substantial tuning of the gauge coupling. 
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1 Introduction 

Cosmic Inflation is considered by most theoretical cosmologists as an integral part of the 
so-called concordance model of cosmology, which complements the standard hot big bang. 
Inflation not only overcomes the flatness and horizon problems but, crucially, it provides 
an elegant mechanism for the generation of the primordial curvature perturbation, which 
is responsible for the origin of structure in the Universe [1]. Because of this, in the last 
25 years there has been a massive effort to design compelling models of inflation based on 
realistic particle theory. This endeavour was accompanied with the influx of observational 
data on cosmological perturbations with ever increasing precision, which have the power to 
discriminate between inflation models. However, despite the fact that observations are now 
putting substantial pressure on inflation model-building, we are still far from ascertaining 
which inflation model is the most likely, since no candidate has been found with a clear 
advantage over others. 

In the early years of inflation model-building emphasis was put on economy. This gave 
birth to the inflationary paradigm which considered that inflation and the cosmological per- 
turbations are determined by a single degree of freedom, a scalar field called the inflaton. 
This was not only because inflation had to outperform the rival paradigm of cosmic strings 1 
but also because it was the simplest choice, which allowed a handle over calculations and 
provided inflation with predictive power. Conceptual problems regarding initial conditions 
were overcome by the no-hair theorem (which put the issue beyond our reach) and by the 

1 Cosmic strings had an important advantage over inflation in those days, because they could explain the 
magnitude of the cosmological perturbations ~ 10 -5 . 
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idea of eternal inflation [2], which is an initial phase allowed in many inflationary models. 
However, insisting on a single degree of freedom resulted unavoidably to massive fine-tuning, 
which was still far less severe than that of the horizon and flatness problems, but nevertheless, 
it was found to plague the vast majority of inflation models. 

One of the sources of fine-tuning was the apparent requirement of super-Planckian val- 
ues of the inflaton field in large-field inflation models. 2 Super-Planckian field displacements 
implied that non-renormalisable terms in the scalar potential could not be ignored. In fact, 
there was real danger of them blowing up, rendering the theory out of control. 

A seminal step to overcome the fine tunings of single-field inflation was made by Linde 
in ref. [3], where he introduced a second scalar field (the so-called waterfall field) interacting 
with the inflaton and causing the termination of the inflationary phase. The model was named 
Hybrid Inflation and is one of the most promising types of inflation models that remain into 
effect today. In Hybrid Inflation the inflaton is relieved from one of its responsibilities (to end 
inflation) and this is why the tuning requirements of the model are reduced, at the expense 
of adding another degree of freedom. This enabled the inflaton field to stay sub-Planckian 
during inflation. Moreover, since inflation had to take place at energies near the breaking of 
grand unification (to satisfy the observational requirement on the magnitude of the curvature 
perturbation), the waterfall field could be identified with the Higgs field of a Grand Unified 
Theory (GUT). 

However, the original Hybrid Inflation model suffered also from the problem of exces- 
sive radiative corrections, which lifted the flatness of the inflaton direction. To harness this 
effect an additional symmetry was needed. Global supersymmetry (SUSY) was proposed in 
ref. [4]. SUSY offers a multitude of flat directions whose flatness is guaranteed by the non- 
renormalisation theorem. This is why SUSY has been extensively used in inflation model- 
building. In SUSY Hybrid Inflation the radiative corrections are actually beneficial because 
they provide a gentle logarithmic slope for the inflaton to slow-roll. 

However, local supersymmetry, called supergravity (SUGRA), introduced another set of 
fine-tunings because it turns out that Kahler corrections to the scalar potential generically 
give rise to masses of order the Hubble scale H to all scalar fields in the theory [5]. This is 
the basis for the infamous r/-problem of SUGRA inflation, which has to do with the fact that 
the slow-roll parameter rj was pushed to order unity by SUGRA corrections, where 

2 V" 1 /m\2 

with m ~ H being the mass of the scalar field, the prime denotes derivative of the scalar 
potential V with respect to the scalar field and we used that the Friedmann equation during 
inflation is ?>m 2 p H 2 ~ V, since the Universe is dominated by the potential density of the 
inflaton. Note here that the above argument is true for any light scalar field during inflation; 
not only the inflaton. This is a problem because, if a scalar field contributes to the curvature 
perturbation then the spectral index n a receives from 7] a contribution of the form 

5(n s - 1) = 2n . 

Thus, to obtain the approximate scale invariance revealed by the CMB observations, one 
requires ijCl. A second problem has to do with the requirement that the inflaton slow-rolls, 



2 The small field models suffered from another fine-tuning problem since the hilltop in the scalar potential 
had to be extremely flat. 
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since we need about 60 e-folds of inflation to solve the horizon and flatness problems. However, 
SUGRA corrections lift the flatness of the inflaton direction and destabilise slow-roll. Finally, 
a third problem is that, if the scalar field mass is m > particle production is exponentially 
suppressed, so there is danger of not being able to generate the density perturbations using 
scalar fields. 

For Hybrid Inflation, the above does not affect the waterfall field whose mass is generally 
much larger (except momentarily near the end of inflation), but it could potentially affect the 
inflaton direction. Fortunately, when a minimal Kahler potential is assumed, a cancelation 
protects the model from the excessive supergravity corrections, and allows it to escape the 
//-problem [6]. Yet, the price to pay was that a minimal Kahler potential has to be assumed 
because any generic Kahler potential (or any higher order corrections to the minimal Kahler 
potential) would produce a massive 77-problem and render the model nonviable. 

In the meantime the observation of the CMB acoustic peaks lead to the collapse of 
the cosmic string paradigm and put inflation at the centre-stage as the most compelling 
paradigm for the origin of structure. Gradually, model-builders began to move away from 
single-field inflation models of the old days and started designing more complicated models 
with better connection to particle theory. It was soon realised that adding new degrees of 
freedom reduces the fine-tuning requirements of inflation models. Thus, there was a wave of 
new models involving many scalar fields such as double inflation [7], extended inflation [8] 
and more recently assisted inflation [9] (or N-flation [10]), where a cascade of hundreds of 
identical scalar fields is assumed. Another approach mirrored the idea of Hybrid Inflation in 
that it removed the responsibility from the inflaton to accomplish some of its tasks, assigning 
it to another field. In the curvaton paradigm, for example, the inflaton is no more required 
to generate the curvature perturbation through its quantum fluctuations. Instead this task 
is given to some other so-called curvaton field [11], which is unrelated to the physics of 
inflation. Many curvaton candidates were found in simple extensions of the standard model 
of particle physics, while it was demonstrated that inflation model-building was liberated by 
the introduction of a curvaton field [12]. However, being a scalar field, the curvaton also 
suffers from the 77-problem and the same is true for all other proposals which include many 
scalar fields. 3 

Recently a surprising solution to the //-problem of SUGRA inflation was discovered. In 
ref. [16] it was shown that an interaction between the inflaton and a vector boson field affects 
profoundly the inflaton's evolution. Indeed, a new inflationary attractor solution was found 
(called the vector scaling solution) where the vector field backreaction Ba impedes the infla- 
tion's variation in that it reduces the effective potential slope as experienced by the inflaton 
field: \V'g\ < \V'\, where V^' ff = V + Ba and the prime denotes derivative with respect to 
the inflaton field. This allows steep inflation and can overcome the 77-problem by enabling 
slow-roll to take place even if the scalar potential is substantially curved. Furthermore, the 
vector backreaction affects the inflaton equation of motion such that it allows the inflaton to 
undergo particle production even with r\ = 0(1). 

3 The exception is pseudo-Nambu Goldstone bosons (PNGBs), whose flatness is preserved by the remnant of 
the shift-symmetry. Natural inflation [13] exploits this advantage, but unfortunately, to produce the observed 
curvature perturbation, the order parameter of the PNGB has to be super-Planckian, in contrast to quantum 
gravity / string theory considerations. Recently, a promising way out was found in ref. [14] by gauging the 
PNGB inflaton. In this way, vector field backreaction from the axionic coupling facilitates steep inflation, 
which allows the order parameter (decay constant) of the inflaton to be sub-Planckian. This proposal has 
parallels to our work (in the sense that vector backreaction allows steep inflation), even though PNGBs do 
not suffer from the ^-problem per se. For a PNGB curvaton see ref. [15]. 
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In this paper, we demonstrate that the mechanism of vector backreaction not only allows 
long-lasting slow-roll inflation with a steep and curved potential but also protects the spectral 
index of the inflaton's perturbations against excessive contributions from a large rj parameter. 
Thereby, vector backreaction eliminates all aspects of the ^-problem, the excessive tilt of the 
spectrum and the destabilisation of slow-roll for the inflaton. We apply these findings to the 
standard SUGRA Hybrid Inflation model with a generalised Kahler potential, which goes 
beyond the minimal form. We find that even with a non-minimal Kahler potential, Hybrid 
Inflation can be long-lasting and produce a weakly red spectrum of curvature perturbations, 
in agreement with observations. 

Since the pioneering work in ref. [17], which introduced the vector curvaton paradigm, 
cosmic vector fields are increasingly being considered to affect the dynamics of inflation and 
the generation of the curvature perturbation £. Indeed, it was found that, the contribution of 
vector fields can give rise to statistical anisotropy in Q [18, 19]. This is a new observable, which 
cannot be produced by scalar fields alone as the latter cannot generate a preferred direction on 
the microwave sky. In fact, a preferred direction might be hinted by the unlikely correlation 
of the low multiples of the CMB [20]. The observations still allow as much as 30% statistical 
anisotropy in the spectrum of the curvature perturbation [21]. As for the bispectrum, this can 
be predominantly anisotropic even if statistical anisotropy in the spectrum is small [22, 23]. 
The preferred direction of the statistical anisotropy in the spectrum and bispectrum of the 
curvature perturbation is the same, which is a smoking gun for the contribution of a vector 
field to C [23]. 

In ref. [16] it was demonstrated that, when the inflaton affects the kinetic function of 
an Abelian vector boson field, the new attractor solution may render the Universe expansion 
mildly anisotropic (see also ref. [24]) and is such that it generates scale-invariant spectra for 
the vector field components. In refs. [25, 26] it was shown that, anisotropic inflation gives rise 
to statistical anisotropy in the perturbations of the inflaton field which, in turn, pass it over 
to Alternatively, one can employ the vector field perturbations themselves to generate a 
statistically anisotropic Q [18], e.g. if the vector field acts as a vector curvaton as shown in 
refs. [27, 28] (see also refs. [18, 29]). Both contributions are complementary. 

In this paper we apply the above in the generation of statistical anisotropy in the spec- 
trum and bispectrum of the curvature perturbation in SUGRA Hybrid Inflation, by coupling 
the vector field to the waterfall field so that it can play the role of vector curvaton. Such 
a coupling was investigated in ref. [19] to generate statistical anisotropy in £, through the 
end-of-inflation mechanism, since the vector field perturbs the moment when inflation is ter- 
minated by the breaking of grand unification. Our investigation is complementary to that 
work as we study the additional statistical anisotropy that can be generated through the 
vector curvaton mechanism. 

In our paper we consider natural units, where c = h = kp = 1 and Newton's gravitational 
constant is 8ttG = rrip 2 , with mp = 2.44 x 10 18 GeV being the reduced Planck mass. 

2 Vector Scaling Slow-Roll Inflation 

Consider the U{1) gauge symmetry $ — > $e l ^ x \ where the group parameter A(x) is a func- 
tion of spacetime coordinates x, acting on the complex scalar field <I> with unit charge. To 
achieve invariance of the kinetic term |L>^$| 2 , we need to include the Abelian gauge field 
with transformation h^B^ — > HqB^ + d^A(x) via the covariant derivative = du — ih^B^, 
where /iq is the gauge coupling constant. Switching to polar coordinates $ = cf)e t9 /^/2 and 



- 4 - 



defining the (gauge invariant) combination HqA^ = HqB^ — 8^9, we can obtain the following 
Lagrangian density for an Abelian Higgs model 

C = -lf(a)F^F^ + \ (drf) (d^) + \h\^A^ - Vl(</>), (2.1) 

where the field strength tensor is = d^B u — dvB^ = d^A u — d u A^ and V\((j)) is the scalar 
potential for the Higgs field 4>. In this model the gauge kinetic function / is modulated by 
another scalar field a which will drive a period of inflation (this is natural in supergravity 
theories where / is a holomorphic function of the scalar fields of the model). Defining the 
gauge kinetic function as follows 

where we assume that /i(o"o) = ^0 s ° that /(do) = 1 an d the vector field becomes canonically 
normalised when a assumes its vacuum expectation value (VEV) <jq. 

The mass of the gauge field, m = ho4>, is given through the Higgs mechanism by the 
non-zero VEV (fro of the Higgs field (f> which spontaneously breaks the gauge symmetry. 
The spatial components of the physical vector field are given by W% = \/~fAi/a, with mass 
Ma = m/VT = h(a)<j). Hence, once the gauge symmetry is spontaneously broken the physical 
vector field obtains a mass rriA = Ma^Oi^o) = ^O^Oj where we assumed that a — > <tq after 
the phase transition, i.e. the phase transition terminates inflation, as is usually the case in 
Hybrid Inflation. 

It was shown in ref. [17] that, as inflation homogenises the vector field diA^ = 0, the tem- 
poral component vanishes Aq = (if m = 0, we can set Aq = by a gauge choice). Thus, with- 
out loss of generality, the spatial components of the vector field can be lined up in the z axis 
A^ = (0, 0, 0, A z (t)). We can then assume a Bianchi-I background with residual isotropy in the 
plane perpendicular to the vector expectation value ds 2 = dt 2 — a 2 (t)(dx 2 + dy 2 ) — a 2 (t)dz 2 , 
where are the scale factors related to the different spatial directions. We can then define 
the average scale factor a = (afc^) 1 / 3 and the average Hubble rate as H = d/a, where a dot 
denotes a derivative with respect to cosmic time t. The anisotropic stress £ induced by the 
vector field is then given by 

SsE-y-ln/^y (2.3) 
3Hdt \a 2 J 

The coupling of the vector field to the scalar field, through the kinetic function f(cr), 
induces a source term Ba = — a 2 2 f'( a )A 2 /'2 in the scalar field equation 

a + 3Ha + V'(a)+B A ((T,A z ) = 0, (2.4) 

where a prime denotes a derivative with respect to the field a and V(a) is the scalar potential 
associated with the field a. 

As shown in ref. [16], for half of the model parameter space, and for sufficiently flat 
potentials, the system may evolve to the Standard Slow-Roll (SSR) inflationary attractor 
with a as the inflaton. On the SSR attractor the vector backreaction Ba and the anisotropic 
stress S vanish. The vector field energy density pa also vanishes on the SSR attractor and 
hence cannot influence the expansion of the Universe. The SSR attractor provides the typical 
example to the cosmic no-hair theorem [30]. 

For the other half of parameter space, the source term Ba can backreact on the dynamics 
of the inflaton field and lead to novel attractor solutions. In this case the vector backreaction 
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grows and dynamical analysis shows that solutions tend to the Vector Scaling Slow-Roll 
(VSSR) inflationary attractor [16]. 4 On the VSSR attractor the vector field has a non- 
negligible effect on the expansion of the Universe through its non- vanishing energy density and 
hence non- vanishing anisotropic stress. The VSSR attractor, being an anisotropic inflationary 
solution, nicely provides a counter example to the cosmic no-hair theorem [30]. 

On the VSSR attractor, as was shown in ref. [16], the gauge kinetic function scales as 
/(cr) = / att oc a -4 (we notice that because / oc 1/h 2 and / — > 1 at the end of inflation, the 
gauge field remains weakly coupled throughout inflation) . This scaling of the kinetic function 
leads to scale-invariant transverse spectra 5 Vl,R of vector field perturbations [27]. If the 
vector field is massless the longitudinal component and its perturbations are decoupled from 
the theory. In this case, / oc a -4 is all we need to consider, if we are interested in a significant 
effect of the vector field on cosmological scales (such as a contribution to the curvature 
perturbation of the Universe £). 6 Because we are left with the transverse components only, 
particle production of the vector field is highly anisotropic. Therefore, the vector field can 
only contribute subdominantly to £. However for a massive vector field, we have to consider 
also the longitudinal component, which turns out to generate the dominant contribution. To 
obtain a scale-invariant longitudinal spectrum of perturbations V\\ , we also require that m oc a 
[27]. In this case the phenomenology is richer as the vector field spectra depend on the mass 
of the physical vector field at the end of inflation. If the field remains light until the end 
of inflation we have highly anisotropic vector field spectra and the field can only contribute 
subdominantly to £, but it can still generate substantial statistical anisotropy and anisotropic 
non-Gaussianity in the primordial curvature perturbation [27]. If the vector field becomes 
heavy by the end of inflation then its spectra become isotropic and the field can then alone 
generate £ [27]. 

In ref. [16], the VSSR attractor solution, including the energy densities and anisotropic 
stress, was obtained in terms of the dimensionless model parameters defined by 

(2.5) 

4>=o 

where we are considering a massless gauge field, i.e. unbroken gauge symmetry where = 0. 
The backreaction criteria for the VSSR attractor. given in ref. [16], arc simply 

(i |r/|»i, 

Conditions I II |r/| > |A |, (2.6) 

[m A r / >6. 

The VSSR attractor is a stationary point in phase-space if the dimensionless model 
parameters are exactly constant i.e. for exponential forms of the kinetic function and scalar 
potential. However, for more general functional forms of the kinetic function and potential, 
the dimensionless model parameters are time-dependent. Therefore the VSSR attractor is 
really a moving point in phase-space. As shown in ref. [16], additional conditions arise to 

4 A similar stability analysis for the massless version of our model in ref. [16] assuming an exponential 
functional dependence for /(cr) and V(a) was performed in refs. [31] and [24]. 

5 Where 'L' and 'R' denote the left and right transverse polarizations respectively and '||' denotes the 
longitudinal polarization. As this is a parity conserving theory the left and right components are identical. 

6 For a review of the vector curvaton paradigm and statistical anisotropy see ref. [32]. 
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guarantee that solutions are able to reach the moving attractor and be dragged along with 
it. These conditions are given by 

A,B,C,D<1, (2.7) 
where the functions A, B, C and D are defined by 



A = 4 



B 



C 



D 



mp 



- 3 m P 



mp 



- 3 m P 



I", 1 An 



_J_ 

r 2 

An 



+ 



3Tj 



r 2 



r f , A' A(,r / + A r' 



r 2 



+ 



v 

r 2 



2Ai 


Ei 


Ir ? + 


l r 5 




(2- 



(2.9) 



(2.10) 



(2.11) 



2.1 General properties of the VSSR attractor 

In this section we summarise briefly the results obtained in ref. [16] for the VSSR inflationary 
attractor. 

The vector backreaction Ba = —o^ f'{p)A z l1 is not dynamic with respect to the scalar 
field since it is not a function of a. Thus the backreaction may be interpreted as to only modify 
the effective slope of the potential V^ s =V' + Ba- As we are considering gauge symmetry the 
gauge kinetic function / oc l//i 2 , and we assume canonical normalisation / — > 1 at the end of 
inflation, we require that the kinetic function is always decreasing in time, f(t) < 0, so that 
the gauge field remains weakly coupled. Because of this we notice that the backreaction Ba 
always has an opposite sign to the potential slope V'(o). This is because the scalar field rolls 
down its potential, so if V'{a) > 0, then a < and therefore f'(cr) > 0; however if V'(cr) < 0, 
then a > and therefore f'(cr) < 0. This ensures that Ba always has an opposite sign to 
the potential slope V'(a). We conclude that the backreaction (if non-negligible) will always 
reduce the effective potential slope experienced by the inflaton and slow down the scalar field 
as it rolls down its potential. And for the same reasons as described above, the dimensionless 
model parameters An and Tf always have the same sign and thus the combination Anl^ is 
always positive. 

Once the VSSR attractor is obtained, which is a fixed point in phase-space (it may be 
a slowly moving point depending on the dimensionless model parameters), the backreaction 
becomes proportional to the potential slope V'(a). Therefore the effective scalar potential 
slope V^ s as seen by the scalar field is given by [16] 

V + B A ^-r^V. (2.12) 



Considering condition III in eq. (2.6), we observe that the effective potential slope seen by 
the inflaton is reduced. This effect can be used to obtain slow-roll inflation in models with 
potentials that would normally be too steep to support a slow-roll regime. Indeed this at- 
tractor can be used to solve the infamous yy-problem of supergravity (see section 5). On the 
attractor the slow- roll parameters en = —H/H 2 and t/h = —H/2HH are given by 



2Ao 
17 



2Ao 
17 



/ 



(2.13) 
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Therefore a stage of vector scaling slow-roll inflation with £h,??h <C 1 is possible. Hence we 
may consider the slow-roll equations 



3m 2 P H 2 ~ V(a) 



and 



3H& 



(2.14) 



where H is now the inflationary Hubble scale on the VSSR attractor. 

The number of e-foldings of expansion 2V a tt during the period of VSSR inflation, where 
/(o") = /att(o-) oc a -4 , is given by 



1 fjOj) 

iVatt = 7 In 77 r ; 

4 /(CTend) 



(2.15) 



where we denote the field values at the start and end of the attractor as <7j and cr en d respec- 
t ivcly. 

The vector-to-scalar field energy density ratio 1Z does not vanish on the attractor, as it 
would normally do in the standard slow-roll case, but acquires a small non-zero value given 
by [16] 



n _PA^ A IV - 6 



Pa 



r 2 
1 / 



(2.16) 



This non- vanishing vector field energy density induces a small anisotropic stress, given by 
E ~ 2T2./3. The VSSR attractor therefore leads to mildly anisotropic inflation which can also 
source statistical anisotropy in the primordial curvature perturbation £ [25, 26]. 

3 The primordial curvature perturbation from VSSR inflation 

We will now calculate, for the first time, inflationary observables that arise from inflaton 
perturbations generated during the VSSR attractor. At horizon exit of the mode k during 
inflation (in particular, the epoch of horizon exit for the pivot scale k = O.OC^Mpc" 1 [1], which 
we denote by the subscript *), the curvature perturbation generated by inflaton perturbations 
is given by 



Pa 



1 



y3/2 



mi 



V' 



off 



1 V l l 2 \T 



f\ 



30\/2ir 



m. 



(3.1) 



where we used the slow-roll equations in eq. (2.14) together with eq. (2.12). Considering 
that the observed curvature perturbation £ is generated by the slowly-rolling inflaton we 
have Qfj ~ £, where the COBE normalisation is £ ~ 4.8 x 10 -5 [33]. We notice that, on the 
VSSR attractor, the curvature perturbation depends on the inflationary scale and the "slope" 
(modulation rate /') of the kinetic function, through Tf, but not on the slope of the potential 

V. 



3.1 The spectral index and its running 

Considering that the entire curvature perturbation is generated by the inflaton, whose spec- 
trum (assuming that it is a light field) is given by V a = (^p) , we find [1] 



C ~ 4vr 2 1 a 



(3.2) 
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During the VSSR attractor we may use the slow-roll equations in eq. (2.14) evaluated at the 
time when cosmological scales exit the horizon, hence 



1 V(a) fXoTf 



247r 2 mf> e{a) 



G 



(3.3) 



where the slow-roll parameters are defined in the usual way e(a) = (mp/2)(V'/V) 2 and 
rj{a) = mp (V"/V). The spectral index is given by 



1 = 



dlnV, 



a 1 dV, 



dink HV C da 



(3.4) 



where we used d ^ ak = jj-jj^ in the last equality. On the VSSR attractor the spectral index 



becomes 



1 ~ 



2r] — 6e — 2mp 



2(A r / )' 



A r7 



e + mp\ - 



2Ao£/ 
3 F f 



(3.5) 
(3.6) 



This result is completely general for the attractor solution, we have not yet specified any 
model. We recover the standard result n s — 1 = 2ij — 6e for V^g = V, i.e. XoFf = 6, hence 

(A rv)' = o. 

We notice that the spectral index becomes independent of the potential curvature en- 
coded in r/. It is therefore easy to obtain a red spectrum independently of the curvature. In 
effect the spectral index now depends on the "curvature" of the kinetic function (/") through 

F' 

.f 

The current CMB observational bounds from WMAP7 on the scalar spectral index are 
0.953 < n s < 0.981 and n s — 1 = —0.033 ± 0.014 (at la) [33], i.e. they favour a red spectrum, 
although exact scale-invariance is not ruled out. 

The running of the spectral index, given by n' s = on the VSSR attractor becomes 



n Q ~ 



6 



16ei] - 24e z - 2£ 



2 (AoIV)' 
-- mp ^Fj~ 



/2 (A r f V 

6e - 2r? + 2d -m P y " /; 

V 3 1/ 



, 4 2 Ao 
H — m P —- 
3 p F f 



(A r/) 



„ r f( A o r /)' 



r. 



(3.7) 



where ^ = rap (V'V" /V 2 ^j . The standard result n' s = 16er] — 24e 2 — 2^ is again recovered for 
V'g = V', i.e. XoFf = 6, hence (Aol^)' = 0. The result in eq. (3.7) can be written as 



2c 



X F f 



77 - 2e + 2m P 



r" 



r' 



/ 



(3. 



The current CMB observational bounds, from WMAP7, on the running of the spectral 
index with no gravitational waves are —0.084 < n' s < 0.020 (at 95%cf) [33]. 
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3.2 The tensor spectrum 

We should also consider the primordial tensor perturbation which may be generated during 
inflation. Its spectrum Vh and corresponding spectral index ny are given by [1] 

An important quantity is the tensor-to-scalar ratio defined as r = VhjV^- On the VSSR 
attractor we find the following results 

Hence the 'consistency relation' of standard slow-roll inflation, = —r/8, does not hold 
for VSSR inflation. Therefore, if r and nx are observed and do not satisfy the 'consistency 
relation', the evidence could point towards other inflationary attractors such as the VSSR 
attractor. 

The current CMB observational bound from WMAP7 on gravitational waves with no 
running is r < 0.36 (at 95%cf) [33]. 

4 Statistical anisotropy from the vector curvaton 

Clearly, because of the non-vanishing vector field energy density during VSSR inflation 1Z ^ 0, 
we have to consider the possibility that vector field perturbations contribute to the primordial 
curvature perturbation of the Universe [17, 28]. The vector field must be light during inflation 
for successful particle production. But the vector field cannot dominate the energy density 
as it would cause excessive anisotropic expansion in conflict with observations (notice how 
this condition is easily achieved through the VSSR attractor predicting 1Z <C 1). However, if 
the vector field becomes heavy after inflation it begins to oscillate rapidly and behaves like 
pressureless matter [17]. The vector field can then (nearly) dominate the energy density and 
imprint its spectra of perturbations through the vector curvaton mechanism [17]. 

The simplest possibility is that the vector field remains massless during inflation. Indeed, 
consider first that inflation occurs while eft = 0, therefore the vector field is massless. The 
spectra of perturbations generated during inflation for a massless vector field is given by [27] 

Vl,r = V a = and V\\ = 0. (4.1) 

These scale-invariant results for the transverse components of the vector field spectra are 
obtained only if / oc a? or / oc a~ 4 [27]. As seen, the latter condition arises automatically out 
of the VSSR attractor [16]. 

To act as a vector curvaton the vector field needs to obtain a non-zero mass. This is 
achieved by Higgsing the vector field as shown in eq. (2.1). The Higgs field assumes a non-zero 
expectation value at the phase transition which terminates inflation. This in turn generates 
a mass for the vector field. Once the vector field obtains a mass, by the non-zero VEV of 
the Higgs field 4>q which spontaneously breaks the gauge symmetry, it may undergo coherent 
oscillations. The vector field then has a particle interpretation and its decay rate Ta is given 

T a = ^MA = ^ = H dec , (4.2) 
8vr 8vr v ' 
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where the subscript 'dec' denotes the epoch of vector field decay. Let us define the following 
useful parameter 

H en( i 32\/6vr 2 mp^ 

where eqs. (3.1) and (4.2), together with the assumption that -ff en d — (where the subscript 
'end' marks the end of inflation), were used to obtain the second equality. For the gravita- 
tional effect of the vector curvaton not to be suppressed we require that e < 1. From this 
requirement, an upper bound on the gauge coupling constant is obtained 



< 



/ 32V~6ir 2 mp( 

{ 0o|r7(<7*)| 



1/3 



(4.4) 



From eq. (4.1) we see that the vector field spectra of perturbations are completely 
anisotropic Vl,r 7^ V\\ , and therefore cannot be the dominant contribution to £. The con- 
tribution of vector field perturbations will be in generating statistical anisotropy. Statistical 
anisotropy in the spectrum can be parameterized as 



i + 5 (d-k) 2 + 



(4.5) 



where k = k//c, d is the unit vector in the preferred direction chosen by the homogenised 
vector field, g quantifies the statistical anisotropy in the spectrum and the ellipsis denotes 
higher order terms which are not important if g < 1. 

In ref. [21] it was found that g = 0.29 ± 0.03 at the level of 9a. However, the preferred di- 
rection was too close to the ecliptic plane so the authors suspected some unknown systematic. 
Hence, this number can only be considered as an upper bound \g\ < 0.3. The observations of 
the Planck satellite will decrease this bound to \g\ < 0.02 if it is not observed [34]. 

For the vector curvaton case, in the so-called <5iV-formalism [35], it was found in ref. [18] 

that 

Vw - V + -8 

* = f , v^Wl = TT-r- f> ' (4 ' 6) 

where /? = (N A /N a ) 2 < 1, N a = §^ and N A = |N A | where N\ = Therefore we must 

now consider both the inflaton and the vector curvaton contributions to the primordial cur- 
vature perturbation. In the <5iV-formalism, to first order in perturbations, the total curvature 
perturbation is given by 

( = N a 8a + N A 8W. (4.7) 

The total curvature perturbation can also be calculated as the sum of individual curvature 
perturbations multiplied by an appropriate weighting factor 



C = (l - n A j g + n A ( A , (4.8) 

where 

Cl A = PA+PA and Q = -H S A (4.9) 
P + P Pi 



7 If it is to affect the Universe expansion, the oscillations of the vector curvaton should last at least a Hubble 
time before its decay. 



- 11 - 



where the subscript 'i' labels the different components of the Universe content and p and p 
is the total energy density and pressure respectively. We evaluate the vector field curvature 
perturbation when it decays, i.e. when it is oscillating and the average vector field pressure 
PA = [17]. Therefore we obtain 

ftA = =^=-, (4-10) 
7<.dec + n 

where n = 1 + p a /Pa- The vector field may decay either when the scalar field is still oscil- 
lating i.e. p a = (where n = 1) or during the following stage of radiation domination, i.e. 
Pa = Prad = /°rad/3 (where n = 4/3). From eqs. (4.7) and (4.8), we obtain 



NaSW = n A (A = -RdecU 
N a 5a 1 - Cl A Ca ~ n ( a 



(4.11) 



Then, using that the typical perturbation for the fields is given by 8Wl,r = 5a = (H*/2ir) 
we find the statistical anisotropy 



Na T^dec (A 



(4.12) 



N a n (a 

Let us now calculate the individual curvature perturbations Q = —H(5pi/pi). For the infla- 
ton this becomes Q a ~ as the vector field can only contribute subdominanty to generate 
statistical anisotropy. For an oscillating vector field we have [17] 



Ca 6pA 



3pA 



2b A 
3l 



37TW, 

dec 



(4.13) 



dec 

where W osc is the value of the homogeneous physical vector field at the onset of oscillations 
and 'dec' denotes the time of decay of the vector field. Hence, from eq. (4.12) we obtain 

3mr(W osc 

Perhaps we can also estimate the original amplitude of the oscillations W osc . Let us 
assume that, when inflation ends, we enter into the stage of vector field oscillations very 
quickly. Then, we can consider 

(PA)end ~ {PA)osc and {Pa)end ~ (At)osc => "^-end ~ "^-osc • (4-15) 

We can first obtain (pA)end = (^•/ 9 o-) e nd wnere (Po-)end — U(o" en d). The total energy density 
of the vector field is given by a kinetic term plus a potential term pa = Pkin + V A where 
Pkin = a 2 2 fAz/2 an d Va = a^" 2 m 2 A\j2. Now for a heavy vector field that is oscillating we 
have pkin = V A , where the bar denotes the average, and therefore (pa) osc — a^m^A 2 . Then 
assuming (for sinusoidal oscillations) A 2 ~ A 2 /2, where A is the amplitude of oscillations of 
the comoving vector field. We can then find 

(pa)osc = \hl<gwl c , (4.16) 

where we have now neglected the residual anisotropic stress £ as it quickly vanishes during 
vector field oscillations [16] and we considered that W = A/a after the end of inflation. With 
the assumption in eq. (4.15) we find 

2 2U(cT end )^end f41?1 
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then with eq. (4.14) and ~ H end we obtain 



± JWU ^ (41g) 

6 \3nirm P ( J TZ end 



If the inflaton does not decay too early, both fields will be coherently oscillating and 
1Z = constant. The vector field can then decay into a matter dominated background where 

n = 1 and TZ dec ~ TZ end . (4.19) 

However, for a better chance at generating statistical anisotropy, let us consider the case 
where the inflaton decays rapidly into relativistic particles, and the vector curvaton decays 
later into the radiation dominated Universe, hence 

n = ^ and K dec ~ £- 1/2 ft eild . (4.20) 

With this optimal scenario in eq. (4.20), and using eq. (4.3), the statistical anisotropy 
in eq. (4.18) becomes 

191 ~ V 3m P hoC\T f (* m )Y [ ■ ' 

Therefore for Planck detectable statistical anisotropy \g\ > 0.02 [34], in the best case 
scenario from eq. (4.21), an upper bound on the gauge coupling can be obtained 

3m P QTf(a*)\ 
4.1 Anisotropic non-Gaussianity 

The vector curvaton may also generate non-Gaussianity in the primordial curvature pertur- 
bation. 8 The latest CMB observations provide a hint for the detection of the non-linearity 
parameter /nl which characterizes non- Gaussianity f^ al = 32 ± 21 (at la) [33]. For the 
spectra in eq. (4.1) and using eq. (4.6), from ref. [23] the non-linearity parameter was found 
to be 

where W± = |Wj_|, with Wj_ being the projection of the unit vector W to the plane defined 
by the three ki 2 3 vectors which determine the bispectrum. Therefore the non-Gaussianity is 
itself anisotropic, it has directional dependence (given through W±) related to the preferred 
direction in the spectrum. Further more, the amplitude of non-Gaussianity is correlated to the 
amplitude of statistical anisotropy in the spectrum (given by g) . These effects are signature 
predictions of a vector field contribution to the primordial curvature perturbation. 

Using eq. (4.18), the amplitude of non-Gaussianity generated in the sudden-symmetry 
breaking approximation in eq. (4.15) is given by 



In this paper we discuss only the bispectrum, for the trispectrum see ref. [36]. 
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In the best case scenario with eq. (4.20), the amplitude of non-Gaussianity becomes 



/nl=4<7 2 ^. (4.25) 

"J '^end 

Therefore for Planck detectable non-Gaussianity /nl ^ ®(X)i m the best case scenario, the 
following upper bound on the gauge coupling constant must be satisfied 

710 ~ 3888 ^1^(^)1 fej ' ^ 
5 Vector Scaling SUGRA Hybrid Inflation 

We now embed the above vector curvaton into a well motivated model of SUSY GUT Hybrid 
Inflation [4, 6]. As described in previous sections, through the vector-scalar coupling, we will 
obtain a model of vector scaling Hybrid Inflation. 

Consider a simple SUSY GUT model based on the gauge group G = Gsm x U(1)b-l 
where Gsm is the Standard Model gauge group. This model naturally incorporates [4, 6] the 
standard SUSY realisation of Hybrid Inflation. The gauge group G can arise from concrete 
GUT models with larger gauge symmetry. For example the "semi-shifted" [37] inflation- 
ary scenario, which is based on the extended supersymmetric Pati-Salam model, a U(1)b-l 
remains unbroken during inflation which is then spontaneously broken at the end of inflation. 

Consider a conjugate pair of Gsm singlet left handed superfields <I> and <!> with charges 
+1 and —1 respectively with respect to the U(1)b-l g au g e symmetry which thus break the 
U(1)b-l by their VEVs. Consider also a gauge singlet left handed superfield S which will act 
as our inflaton. The flatness of the inflationary trajectory is normally guaranteed by a U(l) 
R-symmetry on S (which only allows terms linear in S in the superpotential W to all orders) , 
however in our model we will consider a discrete Z n symmetry under which S — > Se 27ri / n 
and W — > We 2m / n . Note, in passing, that global continuous symmetries such as the U(l) 
R-symmetry normally considered can effectively arise [38] from the rich discrete symmetry 
groups encountered in many compactified string theories (see e.g. ref. [39]). Allowing for 
non-renormalisable terms, the most general form of the superpotential relevant to inflation 
and permitted by the symmetries is 

oo 

W = S A klk2 {^) kl (S n ) k2 , (5.1) 

fcl,fc2=0 

where Af tl k 2 are coefficients with varying dimensions (for a similar analysis see ref. [40]). 
Expanding the superpotential we find that for n > 3 

W = kS — M 2 ) + "non-renormalisable terms", (5-2) 

with Aqq = —kM 2 and Aiq = K and the parameters k and M are made positive by field 
redefinitions. We will take the view that terms in the superpotential with n / are mp 
suppressed corrections to the renormalisable superpotential given above, for example the 
coefficient Aqi = —ckM 2 /m p for some expansion coefficient c of order unity. 

Through the remainder of this paper, unless stated otherwise, we will associate the scale 
M ~ MquTi where Mqut = 2.86 x 10 16 GeV is the scale of grand unification, so that the 
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spontaneous symmetry breaking corresponds to the breaking of grand unification. The SUSY 
minimum is at (S) = and (<$) = (<J») = M. 

Models in supergravity are defined by three fundamental functions; the superpotential 
W, the Kahler potential K and the gauge kinetic functions f ab . The scalar potential in 
supergravity has the form 



V = e K/m 



1 \w\ 21 

13 i mj, 



+ lj2^ Re f^T 1 h a h b D a D b (5.3) 



2 

a,b 



where 



d 2 K dW W dK dK 

The subscripts a,b, - ■ ■ label the generators T a of the gauge group with gauge couplings h a . 
The £ Q are Fayet-Iliopoulos D-terms, which can only exist for U(l) gauge groups. In our 
model, only the gauge kinetic function associated with the U(1)b—l, which we simply write 
as /, will not equal to unity. All kinetic functions associated with other gauge symmetries 
will be set equal to one. And for our model we have $j = (S, <]?, $). The D-flatness is given 
by |$| = |$| i.e. $* = Qe td where we choose 9 = so that the SUSY vacua are contained. 
Bringing the fields $ and <£> on the real axis by appropriate U{V)b-L transformations, we 
write <J> = = (j)/2 where (f> is a normalised real scalar field. It will also be useful to define 
another real scalar field a where \a\ = V2\S\ (see below). The supersymmetric minimum is 
then found at a = o"o = and <j) = <j)Q = ±2M. 

The effective mass-squared m\ for the field (f) is given by 9 



m: 



no 



k 2 (\S\ 2 -M 2 ), (5.5) 



hence the field becomes tachyonic at the critical value a = a c = ±\/2M. This is the field value 
at which the waterfall takes effect. For |a| > \u c \ the potential is minimised at $ = $ = 
(i.e. = 0) and inflation is driven by the false vacuum energy density k 2 M 4 . Inflation is 
only stable for values |<r| > |<r c |, we can therefore assume that inflation ends abruptly at the 
waterfall. 

The Kahler potential is a real function of the superfields and must be invariant under 
the symmetries of the theory. In our model it is therefore a real function of \S\ 2 , |<&| 2 , |^| 2 , 
and S n . We expect the Kahler potential to be an rap suppressed expansion of the superfields. 
On the inflationary trajectory where $ = $ = 0, it is clear that the matrix ify* becomes 
diagonal. We also notice that the terms F§ = F| = on the trajectory, therefore the only 
terms in the Kahler potential that will contribute on the trajectory are 



°° IS1 2fcl r i 

K = E L + n k2 -2 ^ (S^ 2 + h - C ■ ( 5 - 6 ) 



^^2=0 m P 



9 If the waterfall field is coupled to the gauge field as well, as we assume in this paper, there is also a 
contribution from the gauge field to its mass-squared. According to eq. (2.1), this contribution is ~ (hoW) 2 , 
where W = |W| with Wi = \fjAi/a and / = 1 at the end of inflation. Assuming that W is no larger than 
the value of the physical vector field condensate after the end of inflation, then pw < pint = k 2 M 4 , suggests 
that W < kM/Hq, where pw ~ (hoMW) 2 is the density of the vector field after the phase transition which 
ends inflation. This condition guarantees that the contribution of the vector field to the tachyonic mass of 
the waterfall is subdominant. However, perturbations of the vector field can perturb the phase transition and 
generate a contribution to £ as shown in ref. [19]. We do not consider this mechanism in this paper. 
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The afc 1 fc 2 are dimensionless coefficients of order one. Hence K = \S\ 2 — (a / 4)\S\ 4 / m 2 p + • • • , 
with aoo = 0, a\Q = 1/2 and 020 = —a/8 where \a\ ~ 1 is a real parameter. 

The scalar potential, given in eq. (5.3), on the inflationary trajectory can then be pa- 
rameterized by 

00 I Q|2fci / qn\k 2 

V = #M* £ Pk J-^U +h.c. (5.7) 

fc 1 ,fe 2 =0 m P 

where the Pk±k 2 are dimensionless coefficients that are functions of the superpotential co- 
efficients Ak x k 2 an( i the Kahler potential coefficients a/ Cl fc 2 . Switching to polar coordinates 
S = |5|e l1? = \a\e % ® /V%, the scalar potential, calculated up to terms suppressed by m p , is 
given by 

( / \ 2 / \ 4 

a 



n ^ = 0)^M<jl + |(^-) +Si (£) +27 („ + l)(^L) ccrf + .-.J, 

(5.8) 

where 8/3 = 1 + 7a/2 + 2a 2 — 18(030 + c.c) and 7 = + aoi — an and we assumed 7 to be 
real. The phase will be homogenised by inflation and will assume a constant value. Minimis- 
ing the potential with respect to the phase we find that for 7 < 0, ni9 = 2nk for integer k. Then 
applying I times the generators of the Z n symmetry we obtain n-d = 2nk + 27r/, and appropri- 
ately choosing / we can set $ = 0. For 7 > 0, the potential is minimised with nt? = (2k + l)ir. 
Then applying I times the generators of the Z n symmetry we obtain wd = (2k + l)n + 2ttI, 
and again appropriately choosing / we can set •& = ir/n. 

Looking at the potential in eq. (5.8) we observe that the extra terms arising from the Z n 
symmetry (as opposed to the continuous R-symmetry normally considered) are substantially 
subdominant for n > 3. This potential reduces to the known potential for the canonical 
Kahler potential case where a = 0. 

The one-loop radiative corrections to the scalar potential on the inflationary trajectory 
which are calculated by using the Coleman- Weinberg formula [41] are given by 

AFl _ taop .(^( 21n ^ + /oW ), „ 9) 

with f c (x) = (x + l) 2 ln(l + 1/x) + (x - l) 2 ln(l - l/x) where x = a 2 /2M 2 and TV is the di- 
mensionality of the representation to which $ and belong (in this model N = 1), and Q is 
some renormalisation scale (we note that the slope will not depend on this scale) . In standard 
SUSY Hybrid Inflation, these radiative corrections are necessary to generate a slope on the 
inflationary valley [4] when considering the minimal Kahler potential. 

The non-minimal Kahler potential above induces non-canonical kinetic terms for the 
scalar fields. However, our dynamical analysis discussed in the previous section which leads 
to the VSSR attractor [16], assumes canonically normalised fields. Therefore we have to 
consider sub-Planckian field values a < mp when the backreaction takes effect so that the 
fields are approximately canonically normalised and we can trust our dynamical analysis. 

In the inflationary valley at = 0, the inflaton a obtains a contribution to its mass 
squared from the non-canonical Kahler potential, given by V"((j) = 0) ~ 3aH 2 (i.e. r\ ~ a). 
Therefore the field would normally be fast-rolling and the valley could not support standard 
slow-roll inflation unless the coefficient a is suppressed. However, considering the VSSR 
attractor, there is the possibility to obtain slow-roll inflation even without fine-tuning the 
non-canonical coefficient a. 
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The dimensionless model parameter in eq. (2.5) for the scalar potential given in eq. (5i 



is 



Ao 



a 



(-) 

\mpj 



+ 4/3 - -a' 



\mpj 



+ 



k fmp\ 

8^ \~) 



(5.10) 



For a <C mp there is a competition between the first and last term in the expression above 
depending on the value of the relevant parameters. For large coupling constant k ~ 0.1 and 
| a | ~ 1 there is the possibility that the radiative correction term comes into play near the end 
of inflation. For this term to be important during inflation fine-tuning of the parameter a is 
required. This is usually considered as to obtain a red spectrum of inflaton perturbations. 
However, in this paper, we will consider the natural value \a\ ~ 1 and show that a red spectrum 
can still be obtained if the cosmological scales exit during the VSSR attractor. 

Let us assume that the first term in eq. (5.10) dominates over the radiative corrections 
until the end of inflation at a = a c . The condition for this to occur is given by 



K < 4-7T 



M 
mp 



(5.11) 



For \a\ = 1 we find the upper bound k < 0.148. We will later verify that this condition is 
satisfied for a specific model example. As long as this condition is satisfied we can approximate 
the dimensionless model parameter in eq. (5.10) as 



An ~ a 



3 / a 
2 \mp 



(5.12) 



We note that in eq. (5.12), if \a\ ~ 1 we find |An| < 1. It will also be useful to find the first 
derivative of the dimensionless model parameter 



lm P \X' \ = \r ] -2e\^\a + 3^- 1 -a^ (£)\ 



(5.13) 



6 An exponential gauge kinetic function 

The form of the gauge kinetic function is constrained by its holomorphicity and the symmetries 
of the model. The Z n symmetry prevents the appearance of terms in / which are linear in 
S, but combinations S n are allowed. Combinations are also allowed by the symmetries, 
but these will not contribute on the inflationary trajectory. Consider a gauge field with an 
exponential gauge kinetic function of the following form 



f(S r ' 



exp 



17 



exp 



fj£LY 

\V2MJ 



And 



(6.1) 



where the vector field becomes canonically normalised with / — > 1 as the inflaton settles into 
the SUSY vacuum an = at the end of inflation. For 7 < 0, where we can set •& = we choose 
q to be a real positive parameter and for 7 > where we can set •& = ir/n we choose q to be 
a real negative parameter so that the exponent is always positive 



/(CT) = e M(MA/2M)"_ 



(6.2) 
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The dimensionless model parameter defined in eq. (2.5) is then given by 



r /(o") = kl 



n 



a/3 fmp\ ( a \ ( \a\ 
2 Km) \\a\J VV2M 



n-l 



(6.3) 



For this model, the dimensionless model parameter is not a constant for n/1. The 
existence condition [16] for the VSSR attractor demands that Tf and Ao have the same sign, 
hence we require that the Kahler potential parameter a > 0. We now choose a > for 
definiteness without loss of generality. 

We must now identify the conditions on the model parameters which will lead to the 
VSSR attractor, namely conditions in eqs. (2.6) and (2.7). For \q\ ~ 1, a > a c and n > 3 we 
find Tf^$>l, hence condition I in eq. (2.6) is readily satisfied. Considering the approximation 
in eq. (5.12) with a ~ 1 we have Ao < 1, therefore condition II in eq. (2.6) is also readily 
satisfied. Then requiring that condition III in eq. (2.6) is satisfied right up until the end of 
inflation at a = <r c we obtain the following condition on the model parameters 



\q\an > 4. 



(6.4) 



In this model the VSSR attractor is a moving point in phase-space (because of Ao = Ao(c) 
and Tf = Tf(a)), therefore we have the additional conditions in eq. (2.7) to ensure that 
solutions are able to reach the moving attractor and be dragged along with it. For this 
model, in which T'f/Tf = (n — l)/cr and A /Ao = 1/c, we find from eqs. (2.8)-(2.11) that the 
two tightest conditions are A, C < 1, i.e. 



A 



4\/-mp 



r' 



r 2 
L f 



< i 



and 



C 



4\/-mp 



r 2 



T f (A0I/-6) 



< 1. 



(6.5) 



For A < 1 to be satisfied up until the end of inflation we require that |g| > 8(n — 1) /3n which 
can be easily satisfied even for large n giving \q\ > 8/3. Then for condition C < 1 to be 
satisfied up until the end of inflation, and assuming condition A < 1 is satisfied, we obtain 
the following bound on the parameters 



Ian — 4 > -an. 

I „ 3 



(6.6) 



The condition C < 1 simply tells us that we cannot be too close to the bifurcation point at 
XoTf = 6 (see ref. [16]) as this is where the VSSR attractor disappears and we enter into a 
stage of fast-roll inflation. The bound in eq. (6.6) can be written as \q\ > 8/3 + A/an which 
for large n goes to the same limit as that of condition A < 1 i.e. \q\ > 8/3, but otherwise it 
is a stronger bound on parameter q. The bound in eq. (6.6) is also stronger than the bound 
in eq. (6.4) and therefore guarantees that inflation ends before condition III is violated. This 
bound should be made strong to guarantee that the system remains at the VSSR attractor 
right up until the end of inflation. For a given q and n, the bound on a becomes 



n 



(6.7) 
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6.1 Properties of vector scaling slow-roll inflation 

We will now assume that the VSSR attractor is obtained. From eq. (2.15), the e-foldings iV a tt 
generated during VSSR attractor is given by 



AU. = ^ 



\y/2M 



1 



mp 



V2M ( AN ; 



mp 



'att 



+ 1 



1/n 



(6.8) 



where <jj denotes the field value at the start of the attractor and the field value at the end 
of the attractor is given by cr en d — o'c- Requiring sub-Planckian initial field values cr, < mp, 
an upper bound on the total number of e-folds is obtained for a given q and n. For \q\ = 3 
and n = 3 we find that N att < 10 5 which is many more e-foldings than the number needed to 
solve the horizon and flatness problems. As an example, evaluating the above for prompt re- 
heating where N att — 67 to solve the cosmological problems, with n = 3 and \q\ = 5 we obtain 
ai/mp ~ 0.06 and therefore we are well into the sub-Planckian regime and the assumption 
that the attractor solution is obtained by this field value is justified. Similarly, the field value 
when cosmological scales (in particular the pivot scale) exit the horizon at AT* ~ 60 is given 

by 

'4AL \ 1/n 

*+l] . (6.9) 



(T, 



V2M \ |9| 

For the gauge kinetic function under consideration, we find from eq. (2.12) the reduced 
effective scalar potential slope V^/V ~ (A/\q\an)(\/2M /a) n <C 1. The slow-roll parameters 
given in eq. (2.13) become 



2a f a \ 2 / V2U\ 

e H ^ -r-j- and ?? H ^ £H + 

\q\n \mp J \ a I 



2(2 - n) ( V2M" 



\q\n 



(6.10) 



hence €h, ?7h *C 1 for n > 3 and mp > a > a c = \/2M. 

The attractor solution for the vector-to-scalar energy density ratio 1Z, given in eq. (2.16), 
for this model becomes 



(>L\*(' 



(qn) 2 \m P J \y/2M 



2(l-n) 



an 



(— Y 

\V2Mj 



(6.11) 



Thus, we see that for n > 3, 1Z is increasing with time. Therefore, we must have 1Z <C 1 at 
the end of inflation so that the vector field remains safely subdominant throughout VSSR 
inflation. The value of 1Z(a) at the end of inflation is 



\mp J 



\q\an 



(qnf 



> 16 / M V a > 64 / M x 
~ 3 Imp / q 2 n ~ 3 \mp 



(qn)' 



(6.12) 

where we used eqs. (6.6) and (6.7). With n = 3, \q\ =3 and a = 4 which saturates the bound 
in eq. (6.7), we find that 7£ e nd — 1 x 10~ 4 . For large n the ratio decreases as 7£ e nd oc 1/n. 

6.2 The primordial curvature perturbation and gravitational waves 

From eqs. (3.1) and (6.9), to leading order in the potential in eq. (5.8) i.e. V(o*) ~ k 2 M 4 , 

(n-l)/n 



we obtain 



K\Q\n 



20VQir \mp J V \Q 



.13) 
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hence the COBE normalisation with M = Mqut> N* = 60, n > 3, and \q\ > 3 leads to the 
upper bound on the coupling constant n < 1.5 x 10~ 3 . This bound for the coupling constant 
k readily satisfies the condition in eq. (5.11) for a = 4 i.e. k < 0.30. Therefore the approxi- 
mation used in eq. (5.12), where the radiative corrections to the potential remain negligible 
until the end of inflation, is well justified. The COBE normalisation is satisfied for exactly 
M = Mqut whereas in the standard SUSY Hybrid Inflation model [4, 6] the maximal M 
which can be achieved is somewhat below the SUSY GUT VEV. This is an additional at- 
tractive feature of our model as the spontaneous symmetry breaking which ends inflation 
corresponds exactly to the breaking of grand unification. 

We observe that the second term in the square brackets in eq. (3.6) for the spectral 
index is positive, therefore a red spectrum of perturbations is obtained. For this model, using 
eqs. (5.12), (6.1) and (6.9), we find the following result for the spectral index 



8 /4iV* 



\q\n 



+ 1 



M 
mp 



4iV* \ 2/n , 



.14) 



For n > 3 and ~ 60 we observe that the spectral index depends very weakly on the model 
parameters q and a for \q\, a ~ 1. In fact the result above is well approximated by the compact 
formula 

2(n - 1) 



n, ~ 1 



.15) 



We find that for n = 3 the spectral index is n s ~ 0.978 and for n> 1 we obtain n s ~ 0.967, 
which fit very well within the (la) CMB observational bounds from WMAP7 on the spectral 
index at 0.953 < n s < 0.981 [33]. For large enough n, we can therefore obtain the observa- 
tionally preferred central value for the spectral index n s = 0.967. 

The running of the spectral index, given in eq. (3.8), for this model becomes 



24 



rj-2e- 2n(n - 1) 



mp 

(7* 



(6.16) 



Using eq. (6.9) and with e ~ (a 2 /2)(ff*/mp) 2 and rj ~ a for a ~ 1 and n > 3 the running is 
well approximated by 

2(n-l) 



n 



(6.17) 



We find that for n = 3 the running of the spectral index is n' s ~ —3.6 x 10 and for n>lwe 
obtain n' s ~ —5.4 x 10 -4 , which satisfies the CMB WMAP7 observational constraints on the 
running of the spectral index with no gravitational waves at —0.084 < n' s < 0.020 (at 95%cf) 
[33]. 

The tensor-to-scalar ratio, given in eq. (3.10) and using eq. (6.9), for this model becomes 



256 
(qn) 2 



M 
mp 



AN \2(l-n)/n 



(6.18) 



With ./V* = 60, n > 3 and \q\ >3we find that the running is r < 1 x 10~ 6 , which satisfies the 
current CMB observational bound from WMAP7 on gravitational waves with no running at 
r < 0.36 (at 95%cf) [33] but is probably too small to ever be observable. 
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6.3 Statistical anisotropy and anisotropic non-Gaussianity 



We will now consider the vector field to act as a vector curvaton. A vector field contribution to 
the curvature perturbation will induce statistical anisotropy and anisotropic non-Gaussianity 
[18, 23]. For the gravitational effect of the vector curvaton not to be suppressed we require 
that £ < 1 where e is defined in eq. (4.3). From eq. (4.4) and using eq. (6.9), the following 
upper bound on the gauge coupling constant may be obtained 



With N* = 60, n > 3 and \q\ > 3 the upper bound on the gauge coupling becomes ho < 0.05. 
The gauge coupling for SUSY GUT is around ho ~ 0.7 and therefore we expect the gravita- 
tional effect of the vector curvaton to be suppressed somewhat. 

The statistical anisotropy induced in the spectrum, given in eq. (4.21) and using eq. (6.9), 
for this model becomes 



where TZ en d is given in eq. (6.12). In eq. (4.22) a bound on the gauge coupling is given for 
Planck [34] detectable statistical anisotropy. In this model the bound becomes ho < 1 x 10~ 4 
for parameter values a ~ 4 with TV* = 60, n > 3 and \q\ > 3. This coupling is far smaller than 
the SUSY GUT gauge coupling at ho ~ 0.7 and therefore observable statistical anisotropy 
through the vector curvaton mechanism is not possible in this model unless the gauge coupling 
is highly fine-tuned. 

The amplitude of non-Gaussianity is given by (see eq. (4.25)) 



where g is given above, e in eq. (4.3) and lZ en d in eq. (6.12). For detectable non-Gaussianity 
/nl ^ C(i) we have an additional bound on the gauge coupling given in eq. (4.26). In this 
model the bound becomes ho < 3 x 10 -10 for parameter values a ~ 4 with A* = 60, n > 3 and 
\q\ > 3. Clearly this gauge coupling is far too small and therefore observable non-Gaussianity 
form the vector curvaton mechanism is very unlikely in this model. 10 

7 Conclusions 

We have demonstrated that, in an inflationary model, when the inflaton also modulates the 
kinetic function of a vector boson field, the backreaction to the inflaton's variation is such 
that allows steep inflation despite sizable Kahler corrections to the scalar potential and also 
produces a mildly red spectral index of inflaton perturbations. In that respect it eliminates 
the 77-problem of inflation in supergravity (SUGRA). 

We have applied the above to a model of SUGRA Hybrid Inflation, where the waterfall 
field is taken to be the Higgs field of a Grand Unified Theory (GUT). The vector field is 
taken to be one of the supermassive GUT bosons, which becomes massive at the GUT phase 

10 However, the contribution to statistical anisotropy in the spectrum and bispectrum from the end-of- 
inflation mechanism may be more substantial [19]. 
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transition that terminates inflation. The gauge kinetic function of this field depends on the 
inflaton and is therefore modulated during inflation. We have shown that slow-roll inflation 
can take place with a generic Kahler potential, that includes higher order corrections beyond 
the minimal form, despite the fact that rj = 0(1). Moreover, we have shown that a red 
spectrum of perturbations is attained, in excellent agreement with observations. Indeed, 
assuming an exponential gauge kinetic function, we have shown that one can obtain for the 
spectral index n s ~ 0.97 — 0.98 with rj ~ 1, with negligible running and tensor fraction. 

It is interesting to note that, in SUSY models based on G*sm x U(1)b~l such as ours, 
the baryon asymmetry of the Universe can be easily generated via a primordial non-thermal 
leptogenesis [42], which takes place during the direct inflaton decay to light particles as shown 
in ref. [43]. The cold dark matter in the Universe can, in principle, consist of the lightest 
neutralino as in many SUSY theories. 

In principle the vector field can also contribute to the curvature perturbation £ either 
through the end-of-inflation mechanism [19] or directly, if it acts as a vector curvaton after the 
end of inflation [17, 32]. We have looked into the latter possibility but we have found that the 
contribution to statistical anisotropy in £ becomes important only when the gauge coupling 
is fine-tuned to unnaturally small values. This is because the attractor solution reached by 
the system during inflation is such that the vector field contribution to the energy density is 
rather small. As a result, a long period of vector field oscillations is required after the end of 
inflation for the vector field to become a significant fraction of the density budget. This can 
be achieved only at the expense of a small gauge coupling. For realistic values of the gauge 
coupling the decay rate of the vector field is not very small which does not allow for a large 
period of vector field oscillations after inflation. However, this does not mean necessarily that 
the vector curvaton mechanism is hopeless in generating observable statistical anisotropy in 
this model. Indeed, the above problem may be overcome if the magnitude of the vector field 
condensate is increased at the end of inflation due to parametric resonance effects. We are 
investigating this possibility in ref. [44] and so far the results are promising. 

The vector field can generate statistical anisotropy through other mechanisms as well. 
For example, statistical anisotropy in £ can be also due to the mild anisotropisation of the 
Universe expansion [25, 26], also a feature of the vector scaling solution [16]. The end-of- 
inflation mechanism is a potentially more efficient way to introduce statistical anisotropy in ( 
due to the vector field perturbations. Since the vector field is coupled to the GUT Higgs field, 
its perturbations modulate the effective mass of this Higgs field. This, in turn, modulates the 
critical value of the inflaton, which triggers the GUT phase transition and terminates inflation 
[19]. However, the vector field contribution to the mass of the Higgs field depends on the zero- 
mode value of the vector field during inflation, i.e. before the phase transition. It is not clear 
how to evaluate this because, during inflation, the vector field is massless (the GUT symmetry 
is restored) so the theory is gauge invariant. This means that one can always perform a gauge 
transformation of the form A — > A + C, where C is a constant spatial vector of arbitrary 
magnitude. Thus, the value A of the zero- mode is gauge-dependent, i.e. it is not simply an 
environmental quantity as is the expectation value of the scalar curvaton during inflation. 
Finally, another method to introduce statistical anisotropy in £ by the vector field is through 
inhomogeneous reheating, since the effective mass of the GUT Higgs (which is modulated by 
the vector field) also determines its decay rate. For the varying kinetic function model, this 
method was touched upon in ref. [28], where it was not found promising. 

All in all we have shown that all the pathologies associated with the infamous //-problem 
of inflation in supergravity are eliminated when considering an interaction between the in- 
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flaton and a gauge field of the theory. Since in supergravity the gauge kinetic function is a 
holomorphic function of the scalar fields of the theory, it is quite natural to expect its modula- 
tion due to the variation of the inflaton. Thus, the appearance of our so-called vector scaling 
solution, is a rather generic phenomenon in supergravity theories with a gauge field content. 
As we have shown, this scaling solution not only allows steep inflation but also generates a 
weakly red spectrum of curvature perturbations in agreement with observations. 
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